  function ssTable=writeSSTable(funcModel,parMat,solveOpt,addSol) 
% function writeSSTable(funcModel,parMat,solveOpt,addSol) 
% Given a matrix of parameters parMat, solve the model and compute the SS
% Will deliver empty output if ssMat is empty 

NC=size(parMat,2); 
for ii=1:NC; 
    [GG,RR,CONS,eu,SDX,ZZ,~,ssMat,~,ssNames]=feval(funcModel,parMat(:, ii),solveOpt,addSol);
    if ~isequal(eu,[1;1])
        display('No Solution')
        if ii==1
            error('First Values does not yield a solution')
        else
            ssTable(:,ii+1)=cellstr(repmat('Nan',size(ssTableLoop,1),1));
            continue
        end
    end
    if isempty(ssMat)==true 
        disp('No steady state to report') 
        ssTable=[];
        return
    end 
    
    if ii==1;
        [NR NS]=size(ssMat);
    end
%     for jj=1:NS;
%         ssTableSub=crtcell(ssMat(:,jj),ssNames,...
%             cellstr(sprintf('Sample %1.0f',jj)));
%         if jj > 1
%             ssTableLoop=merge_cells(ssTableLoop,ssTableSub,1);
%         else
%             ssTableLoop=ssTableSub;
%         end
%     end
%     if ii ==1
%         ssTable=cell( size(ssTableLoop,1) , NC+1 );
%         ssTable(:,[1 2])=ssTableLoop;
%     else
%         ssTable(:,ii+1)=ssTableLoop(:,2);
%     end 
end   
ssTable=[];
% printcell(ssTable); 
